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The paper proposes a self-consistent Green function description of the induced surface supercon- 
ductivity in a disordered three-dimensional topological insulator (TI) coupled to an s-wave super- 
conductor. We recover earlier results regarding the induced spin-triplet p-wave pairing, showing 
that a mixture of p- and s-wave pair correlations appears as a result of broken spin-rotation sym- 
metry on the helical surface of the TI. Unlike the s-wave pairing, the p-wave component is found 
to be suppressed in dirty TIs in which the elastic mean-free path is much smaller than the super- 
conducting coherence length. The suppression is due to the generic nonlocality of the spin-triplet 
correlations, which makes them strongly dependent on the mean-free path in a disordered system. 
In dirty TIs the induced superconductivity is predicted to be predominantly s-wave like. In cleaner 
TIs, however, the p-wave component may reach a magnitude comparable with (but not larger than) 
the s-wave pairing. 
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I. INTRODUCTION 

Topological insulators (TIs) are a novel class of mate- 
rials (see e.g. reviews [l|12) in which the character of elec- 
tron transport varies from insulating in the interior of the 
material to metallic near its surface. Such an atypical 
conduction character originates from specific electronic 
states that, for topological reasons, exist only near the 
surface of the material. The defining property of the TI 
surface states is their spin helicity whereby the spin of the 
charge carrier follows the direction of its momentum vec- 
tor. The spin helicity is preserved in collisions with im- 
purities and, generally, with any nonmagnetic crystal or 
sample defects, due to which TIs have been considered as 
platforms for intriguing applications, ranging from spin- 
tronics to topological quantum information processing. 

One of the recent exciting developments in the TI 
research is the theoretical prediction 3 of unconven- 
tional p-wave superconductivity and Majorana states in 
Tl/superconductor (S) junctions (see also reviews QQj- 
7). In three-dimensional (3D) TIs the superconductivity 
can be induced by depositing an s-wave superconductor 
(e.g. Al, W or Nb) on the surface of the TI material^— 
The unconventional superconductivity arises from the he- 
licity of the surface states, which breaks spin rotation 
symmetry, allowing for mixed singlet s-wave and triplet 
p-wave pair correlations. Theoretical aspects of the su- 
perconducting proximity effect in the TIs have been con- 
sidered in Refs. [I3-HIl30 

This paper addresses the role of elastic impurity scat- 
tering in the superconducting proximity effect on the sur- 
face of a 3D TI. The influence of disorder on supercon- 
ducting properties of TIs still remains largely unexplored. 
On the one hand, the conservation of the spin helicity 
in elastic collisions suggests that the spin structure of 
the pair correlations should remain intact in impure TIs. 
Also, Potter and Leei2 have recently proved the robust- 
ness of the induced superconducting gap in the TI against 



elastic impurity scattering. A related issue has been dis- 
cussed in Ref. [25| for a conventional semiconductor with 
spin-orbit coupling. On the other hand, the Pauli ex- 
clusion makes the spin-triplet pairs inherently nonlocal 
in space and, hence, dependent on a carrier mean-free 
path in a disordered system. One should therefore expect 
suppression of the p-wave correlations in dirty TIs with 
the mean-free path smaller than the superconducting co- 
herence length. This expectation is confirmed below by 
direct calculations employing self-consistent Green func- 
tions of a disordered S/TI bilayer. Concretely, we analyze 
the real-space amplitudes of the induced p- and s-wave 
pair correlations and their evolution from a clean to a 
dirty TI (see also Figs. Q] and [2]). We also show that 
these new results are consistent with the earlier conclu- 
sions^ regarding the robustness of the density of states 




FIG. 1: Amplitude of spin-triplet p-wave correlations F^ p ' 
versus distance R (in units of coherence length of a clean TI, 
£) for different disorder strengths: (a) £/£ — (clean TI), (b) 
t^jl = 5, and (c) £/£ — 30 (dirty TI); £ is elastic mean- free 
path in TI [see also Eqs.l[62ll and (f65l) ]. 
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FIG. 2: Amplitude of spin-singlet s-wave correlations F*-"' 
versus distance R (in units of coherence length of a clean TI, 
£) for different disorder strengths: (a) £/£ = (clean TI), (b) 
I/I = 5, and (c) £/l = 30 (dirty TI) [see also Eqs.ljnU and 

Ell. 



the TI surface is described by a phenomenological sin- 
glet pairing potential. On the other hand, microscopic 
approaches (e.g. McMillan's model 26 ) allow for a more 
general description of the proximity effect in terms of the 
Green functions of the S. We will employ a microscopic 
model close in spirit to McMillan's one^ and its adapta- 
tions to various low-dimensional systems (see e.g., Refs. 
l7lll3il7lil9l f25ll27M 33) ■ Without losing essential physics we 
will treat both the TI surface state and the S film as two- 
dimensional (2D) systems in which electronic states are 
labeled by the in-plane momentum fc. Assuming tun- 
neling coupling between the systems, we can write the 
Hamiltonian of such a bilayer as follows 



# = ~D4^] 



fc.*;' 



H k,k' 



T* 

± k,k' 



T k,k' H k,k' 



A k > 
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(1) 



where Al 



and B k (Ak and Bk) are the creation (de- 



struction) operators of the S and TI, respectively, in the 
Nambu (particle-hole) representation: 



in disordered TIs. 

The subsequent sections give a complete account of the 
theoretical approach adopted in this paper. In Sec. II a 
model for disordered S/TI bilayers is introduced. Section 

III describes a conventional s-wave superconductor which 
is used as the source of the proximity effect. Section 

IV is devoted to induced superconductivity in clean and 
disordered TIs and contains the discussion of the main 
results. 



II. S/TI BILAYER MODEL 

A. Hamiltonian 

We begin by reviewing the superconducting proxim- 
ity effect in a planar interface between a thin singlet s- 
wave superconductor (S) film and the TI surface (see also 
Fig. [3]). Such a hybrid system was first considered by Fu 
and Kane 3 . In their approach the proximity effect on 
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FIG. 3: Schematic of a thin-film superconductor 

(S)/topological insulator (TI) surface contact. 
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The diagonal elements H k k , and H k k , of the matrix in 
Eq. (Q]) are the Hamiltonians of the superconductor and 
the normal TI, respectively, in the absence of the tun- 
neling, whereas the off-diagonal matrix element T kk > de- 
scribes the tunneling coupling between the two systems. 
Below we define explicitly these operators. In the chosen 
basis (|2J| the Hamiltonians H k k , and H k k , are 
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Here h k and h k are the normal-state Hamiltonians of 
the S and N systems with corresponding band structure 
parameters, ^4 S and A N , and Fermi energies, E s and 
E N j^ 3 - A s and \ are the pairing potential and its phase in 
the S, a x , y are Pauli spin matrices (a is the unit matrix), 
and Vk,k' is the disorder potential characterized by the 
correlation function 



(Vk,k'Vk lt k'i) 



Ok—k',-ki+k'u 



(7) 



where a denotes the contact area. Finally, the tunneling 
matrix in the chosen basis is 



Tk ,i 
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(8) 



Since the contact area in such structures is large (of order 
of /um 2 ), we assume that the tunneling coupling randomly 
fluctuates on the surface and treat t ktk > as a random ma- 
trix characterized by the correlation function 



It involves all direct time-ordered products of A k (t) 
(B k (t)) and A\,(t') {B ] k ,{t')) operators, ((...)) denotes 
averaging with the ground-state statistical operator and, 
simultaneously, averaging over the realizations of random 
matrices Vk,k' an d t k ^ k > ©. In Eq. (|T0|) the diagonal 

S N 

entries G and G describe the S and TI, respectively, 

SJV N S 

while the off-diagonal ones G and G are the hybrid 
Green functions due to the tunneling. All entries are 4x4 
matrices in basis ©. 



r 2 

(tk.k'tk 1 .k' 1 ) = Sk-k'.-ki+k'i- 

a 



(9) 



Thus, both tunneling tk,k' and disorder Vk,k' can be 
treated simultaneously and on equal footing, using the 
standard self-consistent Born approximation for the 
Green function of the SN system. 



Using the self-consistent Born approximation with re- 
spect to both Vfcfc' and t kk i, we find the Green func- 
tion Gk,k'(t) — G k (e)S k . k ', where G k {e) obeys the Dyson 
equation 



B. Bilayer Green function 

In order to describe the hybrid S/TI system it is con- 
venient to use a matrix Green function: 



G k (e) = G k (e) + G° k (e)t k (e)G k (e) 



(11) 
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where G k (e) is the Green function of the decoupled 
disorder- free SN, and T, k (e) is the self-energy matrix with 
the following structure in SN space: 



Sfc(e) = 



T z a [u 2 G +T 2 Q \t z o q T z a [u 2 G +1^9 }r z a 
T z a [u 2 G +1^0 \r z a Q T z a a [u 2 G +T 2 G }r z a 

I 



(12) 



In Eq. (p~2j) the terms ex u 2 and T 2 are the disorder and where 
tunneling self-energies, respectively, t z is the Pauli ma- 

• 1 1 1 ^7 S ' N SN,NS 

trix m particle-hole space, and y and y denote 

the momentum integrated Green functions: 



dk s,n 
Gv. . 



(27T) 



j2^ Gk ' (13) 



III. 



MODEL OF THE SUPERCONDUCTOR 



Hu = 



ia y A s e iX 



-ia y A s e %x 



(15) 



For a large Fermi energy E s 3> A s and near the S Fermi 
surface the solution of Eq. (fl4|) is given by (see e.g. Ref. 
0) 



We now demonstrate that Eqs. ([TTj) and dT2^1 recover 
the known results for impure superconductors^ and ex- 
plicitly calculate the Green function Q (TT3l) which will 
be needed later for the analysis of the proximity effect. 
Since the tunneling coupling to the TI has only a minor 
effect on the S, we can calculate the S Green function G k 
from the Dyson equation with T = 0: 



H L 



T zO ]G k 
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. (14) 



G k — 



-ia y A' s e- ix <j (e' s ~ Vs (k)) 



ef-r,Vk)-A' 2 



(16) 



where rj s (k) = hv s (k — k s ) with v s and k s being the 
Fermi velocity and momentum in the S, A' (e,k) and 
e' (e, fe) are the functions of energy e and unit vector 
k in momentum direction on the Fermi surface. These 
functions should be calculated self-consistcntly from the 



equations: 34 
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(17) 
(18) 



where (...) = L -^■{■■■) is the angle averaging over 
the momentum direction k on the Fermi surface and the 
time-scale 



= h/(2TTU 2 V s ) 



(19) 



is related to the disorder strength and coincides with the 
elastic life-time. v s is the normal-state density of states 
(DOS) per spin in the S. 

In what follows we will need the momentum-integrated 
Green function Q (fT5|) which can be obtained from Eqs. 
dHJ), CH]), ([Til) and (UU as 



-mN a 



O og s ( £ ! fc ) i<J y Is ( e > k ) e " 



,(20) 



where the functions g s (e, A;) and / s (e, k) satisfy the equa- 
tions: 



Eqs. (TTTJ) and (|T2|) ]. Up to the T 2 order in tunneling, 

TV 

the equation for G k is 



JV n N 



[e - H k -T T z a G t z <jo - u T z a a Q T z a ]G k = 1,(24) 



where the Green function Q is given by Eqs. (|20[) and 
and H k is the bare surface Hamiltonian: 
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Let us consider first Eq. (|24|) in the absence of disorder: 



[(- - H k -T T z <T g T z a \G k = 1, 



(26) 



where G k denotes the Green function of a clean TI. 
Using Eqs. (f2"U| and (f2"3"| for the momentum-integrated 
Green function Q , we can write Eq. ([26|l as follows 



i/j 



5s(e> fc ) ( A s + ^rfs(^K 



ih 
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(21) 



(22) 



These equations have isotropic solutions g s (e, fe) = g s (e) 
and f s (e,k) — f s (e). In this case Eq. (f2"Tj) becomes 
disorder-independent: g s (e)A s = ef s (e). Solving it to- 
gether with Eq. (|22l) yields the well known result^ 



ffs( e ) 



v/^Aj 



/ s (e) 



A.. 



xA 



A 2 



(23) 



Thus, the momentum-integrated S Green function de- 
pends only on energy through functions 



IV. 



INDUCED SUPERCONDUCTIVITY IN TI 



The induced superconductivity in the TI is described 
by the Dyson equation for the Green function G k [see 



e + tT N (e)-h k -A N (e)ia y e^ 
A„(e)wve-* e + iT N (e) + h"* k 



G k =1, (27) 



where the induced pairing potential A N (e) and the 
spectrum shift iT N (e) are related to the momentum- 
integrated Green functions of the S: 



A N (e)=iT f 3 (e)=iT 

r « ( e ) = T o9s ( e ) = r o 
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(28) 
(29) 
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Here T is the tunneling energy scale. It is the same 
scale that determines the normal-state level broadening 
due to the quasiparticle escape into S with the normal- 
state DOS v s . For large Fermi energy E N ^> A N and 
close to the TI Fermi momentum k N the solution of Eq. 
(|2"T|) is given by 
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where the indices 1, 2 refer to the Nambu space, and func- 
tions e N (e) and rj N (k) are defined by 



e «( e ) = e + * r iv( e ): V N (k) = hv N (k 



(32) 



with v N — A N /H denoting the Fermi velocity in the TI. 
It is also convenient to introduce an alternative form of 
the Green function, G fc , related to (f31"T) by a unitary 
transformation U: 
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sides of Eq. 
— c 2 = 1 and C3 



Up to the cubic terms we find C\ 
3/2, which yields the induced gap 
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(40) 



This equation shows that e g is reduced as the gap in the 
superconductor, A s , becomes smaller. It is also worth 
noting that the energy To can be expressed in terms of ex- 
perimentally accessible parameters of a S/TI interface: 14 
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where the spin and particle-hole sectors are decomposed 
by the direct product (g>, and C k is the Green function in 
the Nambu space only 



Cfc 



e N (e) + 7?„(A0 A„(e)e»* 
A w (e)e-*x e N (e)- ?7jv (fc) 



el(e)- V l(k)-Al(e). 



(35) 



where g„ is the normal-state interface conductance per 
unit area and Uf is the Fermi wave-number of the TI 
surface state determined by the surface carrier density. 
The energy scale Tq can be extracted from the tempera- 
ture dependence of the critical current in short proximity- 
effect junctions^ 



2. Helical s + p - wave pair correlations 



1. Induced gap 

Let us discuss Eq. (|3"Tj) . The single-particle excitations 
in the superconducting TI are described by the Green 

iVO NO 

functions G n k and G 22 k . Both of them involve the 

projector ^ (co + cr ■ k) , indicating that the quasiparticles 

are +1 eigenstates of the helicity <r ■ k, just like in the 
normal state (cf. Ref. 135) . The pole of Eq. (|3Tj) yields 
the equation for the quasiparticle spectrum: 



e 2 N (e)-vl(k)-Al(e)=Q. 
For k close to k N the spectrum is 



£fc ~ ± 



^im 



(36) 



(37) 



It has an induced energy gap e g which satisfies the equa- 
tion e^(e g ) = A^ r (e g ). The latter can be explicitly writ- 
ten as (cf. Ref. [3l|) 



A. 



£o 
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(38) 
where we introduce a dimensionless parameter 7. Assum- 
ing 7 <C 1 for the purpose of this calculation, we search 
for the solution to Eq. (|38| in the form of the expansion: 
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The constants ci, 02,03... are obtained from comparing 
the coefficients at 7,7 2 ,7 3 ... on the left- and right-hand 



I v U i v U 

The off-diagonal entry G 21 k (as well as G 12 k ) in Eq. 
pip is the Green function of the induced superconduct- 
ing condensate. Its structure in spin space reveals the 
symmetry of the induced pair correlations: 



G 2i,fc °c *o-y(^o +cr-k) A N (e) 



(42) 



It is a mixture of the singlet s-wave component (first 
term) and a triplet p-wave component (second term, see 
e.g. Ref. [33). Equation (|4"2")) agrees with the results of 
Ref. [17J for a large Fermi energy. The origin of the mixed 
s- and p-wave superconducting correlations is the broken 
spin-rotation symmetry of the helical surface states. In 
this sense the situation is similar to the mixed singlet- 
triplet intrinsic superconductivity predicted for systems 
without inversion symmetryj 38 i 39 We note that the in- 
duced p-wave component inherits the spin-momentum 
locking cr ■ k of the normal-state carriers. 



B. Disordered TI 

In the presence of disorder it is convenient to recast 
Eq. (|2"41 in the Dyson form: 



Gk ~G k +u T z cr Q g T z a G k , 



(43) 



where G k includes the tunneling and is given by Eq. 
(13T1) . Guided by the decomposition (j33l) and (J34)) we 
seek the solution to Eq. (|4"3"]l in the form 

G N k =UG N k U\ Ol = i(<7o + cr-fe)®G fe) (44) 



where C& is a 2 x 2 Green function in the Nambu space, 
for which we derive the following equation: 
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(45) 



with C fc is given by Eq. (|35|) . The particle helicity results 
in the anisotropic kernel (1 + k ■ q)/2 depending on the 
scattering angle between the directions of the initial k 
and final q momentum states. Since Eq. (1451) contains 
no spin degrees of freedom, it can be solved by analogy 
with the case of the conventional S in Sec. IIII1 which 
yields 
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Here e' N (e, k) and A^ (e, k) are functions of the energy 
e and the momentum direction k satisfying the self- 
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consistency equations [cf. Eqs. (fl7|) and (|18l) ] 



e jv( e ) 



a; ' =a„( C ) 



ifr 


1 + fcq 


e 


TV 


2r N 


2 


v^ 


-A' 2 ' 

JV 


i 

+ - 


h 1 + fe 


g 


A' 

JV 



2r„ 2 ,/e' 2 _ A' 2 

™ V JV JV 



(47) 



(48) 



The bar denotes averaging (...) = Jd*^)...) over the 
momentum direction q of the final state on the Fermi 
surface and the time-scale 



= h/(2TTU 2 P N ) 



(49) 



is determined by the disorder strength u and the normal- 
state DOS per spin in the TI, v N . 

Returning to the full Green function (|44j) we notice 
that it retains the structure of the disorder-free Green 
function (f3"T1) , with A N and e N replaced by A^ ((¥5]) and 



Cfc — 



r N 



12, k 

JV 

22, k 



(cto 



k)(e'+V N (k)) 



{(7q + a ■ k)ia v A' e l 



-i<Jy(a Q + a ■ k) A' e %x -ia y (a + <r ■ k)ia y (e' ' - T] N (k)) _ 



e' 2 - r, 2 (k) - A' 2 

JV 'N v ' JV 



r 



(50) 



1. Solution of the self-consistency equations 

We seek isotropic solutions e' (e, k) = e' (e) and 

A' N (e, k) = A' N (e) for which Eqs. flU} and (gHJ can be 
written as 
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(51) 
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These equations have the same form as in the isotropic 
S [cf. Eqs. ([17} and (JTBJ)] except that the elastic life- 
time t acquires an extra factor of 2 due to the scattering 
anisotropy. Dividing Eq. (jBTj) by Eq. (|52"|) we find 
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(53) 



which upon substituting back to Eqs. (J5TJ) and (I521) yields 
the explicit solutions: 
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These equations along with Eq. (|50"|) provide the self- 
consistent solution for the Green function of a disordered 
TI. 



2. 



DOS 



As an application of the self-consistent Green function 
solution, we now briefly discuss the surface DOS f(e). It 
can be calculated from the formula: 



v{e) 



dk 



r ImTrG llfc (e), 



(56) 



where G n fe (e) is the retarded particle Green function 

which coincides with the G n k (e) block of Eq. (f5"0")) taken 
for e > 0, Im denotes the imaginary part and Tr is the 
trace in spin space. After integration over the momen- 
tum, we have 



v(e) = v N Re 
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where Re denotes the real part. We notice that the 
disorder-dependent factors in the numerator and denom- 
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FIG. 4: Normalized density of states in TI [see Eq. ((58]) ] for 
different interface parameters. 



inator of Eq. ([57]) cancel out, and the DOS assumes 
the usual form for the conventional proximity systems 
(see e.g., Refs. |30| andliol and in Fig. 2]). It features 
two peaks at energies corresponding to the induced gap 
£ g < A s and to the gap of the superconducting material 
A s . The robustness of the DOS with respect to disorder 
is due to the fact that only the isotropic part of the Green 
function ([50]) contributes to Eq. ([56"]) . 



3. Condensate correlation function 

At energies below the induced gap e g , when no single- 
particle excitations exist, the system is a superconducting 
condensate that exhibits pair correlations characterized 
by the anomalous Green function G 21 fe (e) in Eq. (1501) . 
In particular, the anomalous Green function contains in- 
formation on how fast the superconducting correlations 
decay with distance R in real space. The real-space cor- 
relation function can be defined as the Fourier transform: 



G 21 (e,R) = 



dk 



G 2hk (e)e 



ikR 



(59) 



The momentum integration can be done in two steps. 
First, the integration over the absolute value |fc| = k is 
performed. To do so we change the integration variable 
to rj N {k) and extend the lower integration limit to — oo, 
which is permissible for large Fermi energies E N ~^> A' N . 
Then, we do the integration over the momentum direc- 
tion k which is specified by angle <f> in the following equa- 

N 

tions. The resulting expression for G 21 (e, R) is 
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where F^(R) and F^ P \R) are the dimensionless ampli- 
tudes of the s- and p-wave pair correlations given by 



tt/2 

F (s) (i?) = - J cos(k N R 



COS0 X 



x exp 



1 1 

l + 2(i 



Rcoscf) ) dej), (61) 



tt/2 

F {p \R) = - / cos0sin(fc N i?cos(/)) x 



x exp 



1 1 



Rcos(j))d(J3. (62) 



These functions depend on the absolute value of the dis- 
tance R, the Fermi momentum of the TI k N , the coher- 
ence length of a clean TI, £, and the elastic mean-free 
path, £. The latter two are defined by 



s = 



hv» 



V Al(e)-ei(ey 



£ = 



= 2«, 



(63) 



In accord with the Pauli exclusion principle the spin- 
triplet p-wave correlations (1621) are nonlocal, with the am- 
plitude F^(R) vanishing at R = 0. Due to their inher- 
ent nonlocality the p-wave correlations depend stronger 
on the elastic mean-free path £, which leads to the sup- 
pression of the p-wave component in disordered samples. 
This is demonstrated in Fig. [T] for a fixed value of pa- 
rameter k N £ = 10, which correponds to the large- Fermi- 
energy regime. The disorder strength is characterized 
by the ratio £/£ varying from in a clean system with 
£ — > oo (curve a in Fig. [T]) to £/l > 1 in a dirty sample 
with a short mean- free path (curve c in Fig. [I}. The 
comparison of curves a and c in Fig. [T] shows that the 
disorder-induced suppression of the p-wave correlations 
is indeed significant. In contrast, the s-wave amplitude 
([6"T]) is finite for zero distance, F^ s ^(0) = 1 independently 
of disorder (see also Fig. [2]). The effect of a finite mean- 
free path becomes visible only with inceasing R. Still, it 
is much weaker than in the p-wave case. This is clearly 
seen from the asymptotics of Eqs. ([6"T]) and (j6"2"j) at large 
distances R > 2££/(2£ + £): 



F^(R) 



2££ 



nR2£ + £ 



7T R> ? ^ *•' 



F^{R) 



4fc w (2£Q 3 
nR 2 {21 + £f 



*k N e 



TT R 2 ' 

32 h N C 3 

TT R 2 ' 



e«A 



£>1 



(64) 



(65) 



Unlike the s-wave amplitude ([64]) , the p-wave one (|65p is 
proportional to the third power of the mean-free path £ 
in the dirty case £ 3> £. 
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Thus, the main conclusion of this study is that in dirty 
TIs, in which the elastic mean-free path is much smaller 
than the superconducting coherence length (I <C £), the 
induced superconductivity is predominantly s-wave like. 
The p-wave component is greatly suppressed compared 
with the s-wave one. In cleaner TIs the p-wave ampli- 
tude remains always smaller than the s-wave one, albeit 
the difference between them may not be so drastic (see 
curves b for £/£ — 5 in Figs. [T]and[2]). Both compo- 
nents should be observable in the presence of a modest 
amount of impurities. It is worth mentioning that in this 
work the tunneling between a superconductor and a TI 
has been treated as an incoherent process in which the 
electron momentum is not conserved. This case applies 



to spatially nonuniform interfaces. However, similar re- 
sults for the Green functions, the DOS and condensate 
correlation functions are expected for coherent tunnel- 
ing. The self-consistent Green functions derived in this 
paper can be used to study superconducting transport in 
impure TIs. 
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